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In cavity optomechanics, light is used to control mechanical motion. A cen¬ 
tral goal of the held is achieving single-photon strong coupling, which would 
enable the creation of quantum superposition states of motion. Reaching this 
limit requires significant improvements in optomechanical coupling and cavity 
coherence. Here, we introduce a new optomechanical architecture consisting 
of a silicon-nitride membrane coupled to a three-dimensional superconducting 
microwave cavity. Exploiting their large quality-factors, we achieve an optome¬ 
chanical cooperativity of 146,000 and perform sideband cooling of the kilohertz- 
frequency membrane motion to 34 ± 5 microkelvin, the lowest mechanical mode 
temperature reported to date. The achieved cooling is limited only by classical 
noise of the signal generator, and should extend deep into the ground state with 
superconducting filters. Our results suggest that this new realisation of optome¬ 
chanics has the potential to reach the regimes of ultra-large cooperativity and 
single-photon strong coupling, opening up a new generation of experiments. 

In recent years, cavity optomechanics has been used to realise a wide range of exciting 
experiments with mechanical resonators, including achieving exquisite measurement preci¬ 
sion [1-3], strong-coupling between the mechanical and cavity modes [4-6], cooling to the 
quantum ground state [7, 8], for microwave amplification [9, 10], to entangle propagating 
microwave photons with mechanical motion [11], for microwave photon storage [12, 13], to 
observe quantum baekaetion noise [14], to generate squeezed light [15, 16], and to transduce 
photons between the optical and microwave domains [17-19]. These successful experiments, 
performed in the regime of linear optomechanics, were enabled by improvements of the 
coherence of the optomechanical coupling between light and motion. 

In linear optomechanics, the relevant figure of merit for the optomechanical coupling 
is a parameter called cooperativity. Cooperativity combines the optomechanical coupling 
rate g, the cavity decay rate n, and the mechanical decay rate y m , into a dimensionless 
constant C = that quantifies the optomechanical system’s efficiency in exchanging 
photons and phonons [20]. It is similar to the Purcell factor in atomic physics (also often 
referred to as cooperativity) describing the coupling between cavity fields and atoms. An 
advantage of linear optomechanics is that compared to the single-photon coupling rate go, 
the multiphoton coupling rate g is significantly enhanced, given by g = y/N^go where Ad 
is the number of photons used to drive the cavity. The relevant limit for cooperativity in 
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quantum experiments is the so-called quantum coherent limit, in which the cooperativity is 
larger than the number of equilibrium thermal quanta in the mechanical resonator. In the 
quantum coherent limit, an exchange of a photon and phonon occurs faster than the time 
it takes for a phonon to leak into the mechanical resonator ground state from the thermal 
bath. Achieving higher cooperativity that is deeper in the quantum coherent limit in linear 
optomechanics would imply the ability to cool closer to the quantum ground state and the 
preparation of mechanical quantum states with high fidelity. 

Beyond linear optomechanics, the held is striving to reach the limit of single-photon 
strong coupling, in which interaction of light and mechanical motion is coherent at the level 
of a single photon. In this limit, the single-photon coupling rate go exceeds both the cav¬ 
ity decay rate n and the mechanical decay rate y m (go > k, 7 m ). If one could reach this 
limit, one could take advantage of the intrinsic non-linearity of the optomechanical coupling 
at the single-photon level to construct quantum superpositions of mechanical motion using 
optomechanics with classical light. In current experiments, go k, and existing optome¬ 
chanical implementations would require significant improvements in g 0 and/or k to approach 
single-photon strong coupling. 

Here, we present a new optomechanical architecture with large optomechanical cou¬ 
pling that can potentially reach the single-photon strong coupling limit. The design combines 
two highly coherent elements that are applied for the first time in the microwave optome¬ 
chanics domain. The first element is a millimetre-sized, nanometre-thick high-stress silicon 
nitride (SiN x ) membrane, a technology that has demonstrated quality factors up to 1 x 10' 
at cryogenic temperatures [21]. Such membranes have been used extensively in the optical 
domain [14, 16, 22-24] and in transducer applications [18, 19], but have not yet been explored 
as an element in a pure microwave optomechanical system. The second optomechanical ele¬ 
ment is a three-dimensional (3D) microwave cavity, recently popular in the superconducting 
qubit community for their exceptional coherence times [25-27]. By combining these two 
highly coherent elements, we create a new optomechanical platform with large cooperativ¬ 
ity, demonstrate cooling to the lowest mechanical mode temperature reported to date, and 
show that this new system has the potential to scale to couplings significantly beyond the 
state-of-the-art. 

Results 
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Description and characterisation of the device. Fig. 1 illustrates the principle of this 
new 3D optomechanical platform. The mechanical resonator is made from high-stress SiN x 
membrane (Fig. la) that is metallised with an A1 electrode. The cavity itself is an aluminium 
box, Fig. lb, in which electromagnetic fields are confined by superconducting walls in all 
three dimensions. Coupling to the motion of the membrane is achieved by using a flip-chip 
technique to place the membrane on top of antenna electrodes on a separate substrate, 
Fig. lc-e. Fig. If shows an effective lumped-element circuit model of the assembled cavity 
and membrane. Optomechanical coupling results from the modulation of the effective shunt 
capacitance when the membrane displaces, changing the cavity frequency. From finite- 
element simulations, we estimate a single photon coupling rate of g 0 = ( du 0 /dx ) • x zp f = 
27r x 0.36 Hz, where ujq is the cavity mode frequency, x is the mechanical displacement of 
the membrane and x zp { is the amplitude of its zero-point fluctuation. The simplicity of the 
assembly of this 3D cavity architecture also makes it attractive for implementing devices 
such as microwave-to-optical transducers [17-19] by potentially incorporating an optical 
fiber into the superconducting box. 

Measurements are performed in a dilution refrigerator with a base temperature of 
T b = 13 mK (see Supplementary Figure 1 and Supplementary Note 1). Fig. 2a shows a 
measurement of the reflection coefficient |S'n] of the cavity. From a fit to the data, we 
find a total linewidth of k = 2tt x 45.5 kHz, corresponding to a loaded quality factor of 
Ql = 1.1 x 10 5 . From power dependence, we find a maximum intracavity photon occupation 
iV max = 1.3 x 10 8 before the onset of a nonlinear response. The ratio r) = n e /n ~ 0.48 of 
the external decay rate K e and k indicates that the cavity is slightly undercoupled (see 
Supplementary Figure 2 and Supplementary Note 2). 

In Fig. 2b, we characterise the mechanical response of the membrane with the cavity 
using a resonant microwave tone injected at Uq. The thermomechanical motion of the 
membrane generates a peak in the sideband power spectral density (PSD) S(u) of the 
microwave field leaving the cavity at a frequency offset Au = co m from the carrier signal, 
shown in Fig. 2b. We find a mechanical resonance frequency of u m = 2n x 123 kHz, 
consistent with expected fundamental mode frequency of the membrane. A Lorentzian fit 
yields a linewidth of y m = 2n x 3.5 mHz, corresponding to an ultra-high mechanical quality 
factor of Q m = 3.5 x 10 7 , significantly higher than the typical Q m ^ 10 6 of membranes used 
in optomechanical and transducer experiments. 
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Large cooperativity. To quantify the strength of the optomechanical coupling, in Fig. 3 we 
measure the cooperativity C between the mechanical resonator and the cavity. Together with 
initial phonon occupancy of the mechanical resonator nj n , several criteria can be conveniently 
expressed with C, such that for reaching the quantum ground state of motion, C + 1 > n[ n , 
or for reaching the radiation pressure shot noise limit, C > n l m { 1 + (yf-) 2 ). To measure 
C, we use optomechanically-induced transparency (OMIT) [10, 28, 29], which allows one to 
directly determine the cooperativity with no free fit parameters. In OMIT, illustrated in 
Fig. 3a, the cavity is driven by a strong drive tone (oq) while a second weak probe (u; p ) is 
used to measure the cavity response. When driven on the red sideband (tCd = ojq — cc m ), 
a transparency window appears within the broad resonance dip of the cavity reflection 
coefficient [10]. In the limit g <C ft, the lincwidth of the transparency window in |Sii| 2 is 
given by 7 m (C +1) and the peak value by C/{C +1). Fig. 3b shows an example of an OMIT 
measurement with drive-photon number N& = 1.0 x 10 8 : from the broadened linewidth of 
the feature together with its near unity transmission, we extract C = 94,500. Fig. 3c shows 
the extracted C for different drive powers: at the maximum power sustained by the cavity, 
we achieve C max = 1.46 x 10 5 . 


Microkelvin cooling of the membrane resonator. The large cooperativity of our 
optomechanical setup is, in principle, capable of cooling deep into the quantum ground state 
of motion if the mode is thermalised to the temperature of the fridge ( c = 0.015). In 
order to demonstrate the cooling of the resonator, in Fig. 4 we use the spectral density of 


the thermomechanical sideband to directly observe the phonon occupation of the membrane. 
While driving the cavity on the red sideband, the output microwave PSD S(uq) at the cavity 
resonance frequency ccq is given by: 


S(ojq) _ S vv (uj 0 ) 

TiUo QTujJq 




P '^add 


( 1 ) 


where S vv (ujq) is the measured microwave PSD on the spectrum analyser, Q is the net gain 


of the signal path from the output of the cavity to the input of the spectrum analyser, and 


n a dd is the added photon noise quanta from the amplification chain referenced to the output 
of the cavity. By varying the bath temperature T b and measuring S vv (u), we obtain an 
absolute calibration of Q, n a( jd, and n' m (details provided in Supplementary Figure 3 and 
Supplementary Note 3, 4). At the base temperature of the refrigerator, we find that the 
membrane is thermalised to T{ n « 180 mK, corresponding to n| n « 3.06 x 10 4 . From the 
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thermal calibration, we also extract g 0 = 2n x 0.22 Hz, in good agreement with simulations 
(Supplementary Figure 4 and Supplementary Note 5). 

Fig. 4a shows S(uj) for different cooperativities C of the cooling tone. As C increases, 
S(cjq) from the thermomechanical noise peak decreases, indicating that the mode temper¬ 
ature of the mechanical resonator is reduced. At larger cooling powers, however, although 
S(co 0 ) continues to drop, the noise floor of S(co) outside the mechanical bandwidth begins 
to increase, and the peak in the PSD becomes a clip (bottom two panels of Fig. 4a). The 
increase in the noise floor is an indication of noise fluctuations of the cavity field [7]. Due 
to correlations between the fluctuations of the cavity and of the mechanical resonator, the 
spectrum shows a suppression of the total PSD at the cavity frequency [30]. 

In order to extract the mechanical occupation factor in the presence of cavity noise, 
it is no longer sufficient to look only at S^cuo)- In particular, with sufficiently large C 
(and in the absence of other sources of heating that would increase n' m such as losses in 
the superconducting him or in the dielectric membrane or substrate), S^cuo) drops to the 
amplifier noise floor independent of the amount of cavity noise (Eq. 1). This does not, 
however, imply that the mechanical mode is at zero temperature: clue to the hybridisation 
of the mechanical and optical fields, the final mechanical mode occupation in the presence 
of cavity noise is given by (Supplementary Note 6): 


Um ~ (C + l) <+ (C + l) 


n r 


( 2 ) 


where n c is the cavity noise power measured in energy quanta. In order to fold the final 
occupation n m , one must also determine n c . In the limit g <C k, cavity noise appears as an 
increase in the noise floor outside of the mechanical sideband (green dashed line in Fig. 4a) 
with a spectral density (see Supplementary Note 6 for the more general equation of S^cu)): 


S(co 0 + 8u>) 
Ttoj 


= 4?7 ( n c + - ) + n ad d 


( 3 ) 


where (C + l)7 m < du < k. From this expression, one can extract the cavity occupation 
n c and consequently n m for all powers, shown in Fig. 4b. As C is increased, we find that n m 
drops initially to a value of 5.2 ± 0.7, corresponding to a mode temperature of 34 ± 5 /iK, 
beyond which the mechanical occupation begins to increase sharply due to heating from 
cavity noise. The temperature reported here is roughly a factor of two lower than recent 
experiments in the optical domain with silicon nitride membranes [23, 24] due to the very 
low frequency of our membrane. 
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Although the large cooperativity in our experiment should allow us to cool the mem¬ 
brane to n m = 0.2 given the initial thermal occupation, in practice we are limited to 
n m = 5.2 ± 0.7 by cavity noise. In order to cool to lower occupation in future experi¬ 
ments, it is important to identify the source of this cavity noise. The green line in Fig. 4b 
shows the expected n c due to the carrier sideband noise of our microwave signal generator 
(Supplementary Note 7). The good agreement with the observed cavity noise data suggests 
that our final occupation is limited by the spectral purity of the microwave tone used for 
the sideband cooling. 

Discussion 

Having demonstrated the lowest temperature T m = 34 /iK reported to date for a mechanical 
resonator, we analyse the potential of this new implementation to reach the deep quantum 
coherent coupling limit. In the current experiment, the cooling is limited by the classical 
sideband noise of the signal generator. Removing this noise with a tunable superconducting 
cavity with a lincwidth of 10 kHz would already provide sufficient suppression to cool to 
a final occupation of 0.2. A second approach would be to increase go by shrinking the 
capacitor gap: reducing the membrane-antenna gap from 3 /jrn to 100 nm would increase go 
by a factor of 10 3 and C max by 10 6 . A third approach is to improve the cavity linewidth [27], 
also yielding higher cooperativity at lower photon numbers. Finally, combining a smaller 
gap (go ~ 300 Hz) with a better cavity (n ~ 10 Hz) could yield ultra-high cooperativities 
C > 10 12 . Such a fully optimised design would also achieve single-photon strong coupling 
(go > hi, 7m) [31, 32], enabling preparation and detection of non-classical states of motion 
such as Fock states or Schrodinger cat states with optomechanics. 

In conclusion, we have developed a novel optomechanical system coupling the motion 
of a millimetre-sized membrane to a 3D microwave cavity. Exploiting the high coherence of 
the membrane and of the cavity, we achieve a cooperativity of C > 1.4 x 10 5 and perform 
sideband cooling of the millimetre-size membrane to 34 gK, corresponding to a thermal 
occupation n m = 5.2. The scaling of this 3D optomechanical system offers the possibility 
to reach optomechanical couplings far beyond the state-of-the-art, potentially entering the 
single-photon strong coupling regime in which a new generation of quantum experiments 
with mechanical objects would become possible. 
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Methods 

Device preparation. We use commercial SiN x membranes manufactured by Norcada. The 
membranes have the dimensions of 50 nm x 1 mm x 1 mm, and are supplied with a 5 mm x 5 
mm Si frame. We deposit a 20 nm thick him of A1 on top of the membrane without covering 
the clamping edges by using a physical mask. On a separate sapphire substrate, we pattern 
two A1 antenna pads with 80 nm of A1 followed by the deposition of SiN x spacer layer. To 
ensure the membrane does not come into contact with the substrate, a recess of 100 nm is 
etched. The metallised SiN^, membrane is then placed on top of the antenna pads to form a 
capacitor using a vacuum pick-and-place technique. A single drop of 0.1 //I of 2-part epoxy 
is applied on the substrate to attach one corner of the membrane’s Si frame to the substrate. 
Using the depth of focus to locate the vertical position of the membrane and of the bottom 
antenna pads while looking through the membrane with an optical microscope, we estimate 
the gap to be approximately 3 /jm. The gap is much larger than that designed in the spacer, 
most likely due to contamination from dust in the large contact area between the Si frame 
of the membrane and the substrate. The 3D microwave cavity is formed by closing two 
halves of a machined block made out of 6061-aluminium. The inner surface of the cavity is 
polished using a polishing paste, but is not chemically etched. The cavity is assembled by 
screwing the two halves of the cavity together with no sealing mechanism. The dimension of 
the whole cavity is 28 mm(x)x28 mm(y)x8 mm(z), with rounded corners of radius 1 mm 
in the y-z plane. 


Measurement. The bare frequency of the cavity without the antenna substrate and mem- 



brane is 7.4 GHz. The membrane attached to the antenna is placed at the centre of the 
cavity, coupling to the TEno mode. We estimate the total mass of the membrane including 
the A1 layer to be m = 200 ng, corresponding to an estimated quantum zero-point fluctua¬ 
tion of x zp f = 0.6 fm. Measurements are performed in a cryogen free dilution refrigerator. 
To minimise cavity instability during the measurements, the sample is mounted on a mass- 
spring vibration isolation stage on the mixing chamber with resonance frequency of ~ 1 
Hz. During the cooling measurements, the pulse tube is temporarily turned off to minimise 
vibrations of the sample and the cables. Measurements are started 1.5 minutes after the 
switch-off and performed in a ~ 2 minute window before the temperature of the mixing 
chamber begins to increase. 
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FIG. 1. Microwave optomechanics with a 3D superconducting cavity and a millimetre¬ 
sized membrane, a, Electron microscope image of a 50 nm thick SiN x membrane, b, One half 
of the A1 cavity with an SMA connector for reflection measurements. The dimensions of the cavity 
are 28 x 28 x 8 in millimetre, c, Schematic showing the placement of the membrane over the 
antenna pads. The A1 coating of the membrane forms a capacitor with the antennas below, d, A 
complete assembly. The membrane is positioned in the centre of the cavity, supported by a sapphire 
substrate patterned with the A1 antenna, e, An optical microscope image looking from the top 
showing the Al-coated membrane and the underlying antenna pads, f, Effective lumped element 
model of the cavity and membrane. By changing the effective shunt capacitance, the frequency of 
the 3D cavity is modulated by the mechanical motion of the membrane. 


12 
























a 



N 

X 


on 



123.1270 123.1275 

(cl>-6l> 0 )/2tt (kHz) 


FIG. 2. High quality-factors of the microwave cavity and mechanical resonator, a, 

Reflection coefficient |Sii| measurement of the cavity. The loaded quality factor is Ql = 1.1 X 10 5 , 
corresponding to a total decay rate k = 27 r x 45.5 kHz. Internal dissipation rate is kq = 2n x 
23.9 kHz, with Q o = 2 x 10' 5 . b, Power spectral density of the thermomechanical motion of 
the fundamental mode of the membrane uj m ~ 2-7T x 123 kHz with linewdith 7 m = 27r x 3.5 mHz, 
corresponding to Q m = 3.5 x 10 7 . Red: data, Grey: fit. 
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FIG. 3. Large cooperativity measured with optomechanically-induced transparency 

(OMIT), a, Illustration showing the OMIT measurement scheme. The cavity reflection coefficient 
Sii is measured with a weak probe tone (u; p ) while driving the cavity with a second strong tone 
near the red sideband (wd = uq — u m ). A window of transparency appears in the cavity resonance 
(lower panel) with a width set by the mechanical linewdith. b, Extracted cooperativity C vs. 
driving photon number N^. Grey line: expected linear scaling C = 4 ^°^ d . At maximum IVd, the 
cooperativity reaches C max = 1.46 x 10 5 , corresponding to g max = 2ir x 2.5 kHz. 
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FIG. 4. Microkelvin cooling of a millimetre-sized mechanical resonator, a, Measured 
microwave PSD S(u) with red-sideband driving (top to bottom: C = 324, 785, 3810, 94500). 
Black dashed line: noise floor n a dd of the amplifier chain. As the cooperativity C of the cooling 
tone is increased, S{ojq) is reduced and the mechanical linewdith is broadened. Initially, the thermal 
noise of the membrane appears as a peak with decreasing height. At higher C (lower two panels), 
the noise floor outside of the mechanical linewidth (green dashed line) begins to increase above 
7T a dd) indicating the presence of cavity noise. In the lowest panel, the output noise of the cavity is 
suppressed at wo due to correlations of the mechanical and cavity fluctuations. From the thermal 
peak and the cavity noise background level (green dashed line), we extract both n m and n c for all 
cooling powers, shown in b. As a function of C, the mechanical occupancy n m (red circle) drops, 
closely following (blue line) until the onset of cavity noise n c (green squares) limits the final 
occupancy to a minimum of n m = 5.2, corresponding to T m = 34 /iK. The error bars indicate the 
uncertainty of the data points and are calculated with the errors from the Lorentzian fit. The 
green line in b shows the expected n c from the measured carrier sideband noise of the microwave 
generator. 
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Supplementary Figure 1. Experimental setup for reflection measurements. 
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Supplementary Figure 2. Cavity decay rate k as a function of cryostat base tempera¬ 
ture Tb. The line width of the cavity resonance starts to increase above 0.34 K. 


2 





Supplementary Figure 3. Calibration of the initial thermal occupancy by varying 
the refrigerator temperature T\,. The mode temperature of the membrane is thermalised at 
base temperature to ~ 180 mK, corresponding to ~ 3.06 x 10 4 . The error bars indicate the 
uncertainty of the data points, combining the fluctuations due to the temperature and the signal 
level. 
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Supplementary Figure 4. COMSOL model showing a 3D cavity embedding an an¬ 
tenna. 
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Supplementary Note 1: Experimental setup 

The sample is cooled down in a cryogen-free dilution refrigerator to a base temperature of 
12 to 13 mK. Input microwave signals are attenuated at each thermal stage as illustrated 
by Supplementary Figure 1. At base temperature, the signal is injected into the cavity 
via the coupling port of a 20 dB directional coupler. The reflected power collected at the 
output port is amplified by a cold amplifier at the 3 K plate and subsequent amplifiers 
at room temperature and then detected by a network analyser or a spectrum analyser. 
Circulators are placed between the sample and the cold amplifier to isolate the sample from 
noise generated by the amplifier. A cancelling tone is applied using a directional coupler, 
a phase shifter and a continuous attenuator in order to reduce the amplitude of the strong 
drive tone before sending it to the room temperature amplifier. For the optomechanically- 
induced transparency (OMIT) experiment, a network analyser is used to generate the weak 
probe tone ta p and a separate phase locked signal generator provides the strong drive tone tea- 

Supplementary Note 2: Cavity characterisation 

The cavity response |Sii| in Fig. 2a in the main text is fitted with the following equation, 
taking into consideration the finite isolation of the directional coupler [1]: 

ISnMI = ae* + (1 - Q')(l - - - r ) , (1) 

where a is the isolation of the directional coupler. 

Due to the low superconducting transition temperature of aluminium T c ~ 1.2 K, the 
internal cavity quality factor is significantly reduced as the cryostat temperature T h is raised 
above ~0.34 K, and the line width n is broadened, as shown in Supplementary Figure 2. 
This effect needs to be taken into account in the thermal calibration. As a function of 
temperature, we find that n e remains constant. 

Supplementary Note 3: Calibration of the mode temperature 

A thermal calibration is carried out to determine the initial (i.e. without cooling) phonon 
occupancy n' m = ——jd——- of the mechanical resonator, k-g being the Boltzmann constant 
and T' m the initial mode temperature. We send in a carrier signal at ccq to avoid any back- 
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action effects and measure the total power of its mechanical sideband at different cryostat 
temperature Tf,. The total power generated by the thermal motion of the membrane is 
then used to extract the mechanical mode temperature . We use a resolution bandwidth 
(RBW) of 1 Hz. In the limit of RBW> 7 m , the total thermal power can be directly read 
out from the peak height. For T h > 0.34 K, ujq and k changes significantly due to the low 
critical temperature of Al. which has to be taken into account. We plot the converted n' m 
as a function of refrigerator temperature Tf, in Supplementary Figure 3. From the fitting 
line that passes the original point, at the base temperature the membrane is thermalised to 
Tjh ~ 180 mK and n' m « 3.06 x 10 4 . The mode temperature is significantly higher than the 
base temperature of the refrigerator, which is probably an inevitable challenge one is faced 
with when working with a low-frequency (< 1 MHz) mechanical resonator. 


Supplementary Note 4: Power calibration 

The photon number in the cavity corresponding to the input power at the cavity Pj n at 
frequency u is calculated with 


p. 

N = —■ 


K P 


( 2 ) 


hid (f) 2 + (u - u 0 ) 2 ' 

While driving the cavity on resonance, the sideband power due to the thermal motion is 
expressed as 


Pside Pin 


9 o 


k B T 


r Ke \ 2 
< 2 / 


— ^in * g 0 * 2?2 r] 


x zp f ) muj ^ (§) 2 + (v - cj r , b ) 2 (f ) 2 + (a; - a; 0 ) 2 

i 


foe \ 2 
- 2 / 


( 3 ) 


(f ) 2 + (^ — w r>b ) 2 (f ) 2 + (w — u 0 ) 2 ’ 

where go is again the single-photon coupling rate, k B the Boltzmann constant and n m the 
mechanical occupancy. The sideband frequency of interest is either u> r or 07 ,. 

The loss of the input line £ and the gain of the output line Q are calibrated as the 
following: while driving the cavity on resonance (a; = cuo) with power P( n from the signal gen¬ 
erator, the thermomechanical power measured by the spectrum analyser at the mechanical 
sidebands can be written as: 


Pside =P \'n • 9o ■ 2 n m — 




K 


) 2 + 0 ; 2 
V 2 ' 1 m 


■c-g 


( 4 ) 


=K • al • P- 
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The combined value of £ ■ Q is determined by the network analyser. Thermal occupation 
can be measured by sweeping the bath temperature as described in the previous section. All 
the parameters being known in (3, single-photon coupling strength go can be extracted. This 
is equivalent to the single-photon coupling strength calibration the frequency modulation 
(FM) technique described in [8], where the FM peak of the carrier wave provides a side-by- 
side reference of £ • Q. Using this approach, we measure the single-photon coupling strength 
to be go = 0.22 Hz. This further allows us to calibrate the corresponding photon number 
N by using the experimentally obtained photon-enhanced g = goV^N (from cooperativity 
measurements). Subsequently, by rewriting Supplementary Equation 2 as 

kuj (|) 2 + (ca — cuo) 2 ’ U 

we calibrate the total input attenuation £ = 70 dB and correspondingly the output gain 
Q = 73.5 dB. The added noise of the output chain is calibrated to be n a dd = 12 correspond¬ 
ing a noise temperature of approximately 2.6 K. 

Supplementary Note 5: COMSOL RF simulations to estimate go 

To estimate the single-photon coupling rate go of the 3D-cavity-mcmbrane system, we model 
the problem with a “variable capacitor in cavity” structure and numerically calculate the 
cavity frequency u 0 with COMSOL. The following relation is used to derive go'- 

So = || • iC 2 (6) 

where 8 C\ is the change in capacitance by varying the gap between the two antenna rods, 
8 u>o the subsequent frequency shift and 8 C 2 the change in capacitance due to zero-point 
fluctuation x zp f of the membrane. 8 C 2 can be expressed as 8 C 2 = where eo is the 

free space permittivity, A 2 the area of the membrane and d the gap between the membrane 
and the antenna. 

The geometry of the COMSOL model is shown in Supplementary Figure 4. The 
dimension of the 3D cavity is 28 minx28 mmx8 mm. There is a 0.6 mm thick sapphire 
substrate at the centre of the cavity. The antenna is simplified to two A1 rods with a cross 
section of 0.25 mmx0.25 mm, the length of the antenna being 4 mm and their gap is varied 
between 20 and 30 /jrn, changing the parallel-plate capacitance between them. Although 
this geometry is different from the actual geometry of the antenna, it results in comparable 
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frequency pull of c^o, indicating that the capacitance participation ratio in the equivalent 
lumped-element circuit is similar. For a membrane of 1 mmxl mm, a gap of d = 3 /jrn 
and a: zp f = 0.6 fm the resulting coupling strength is go ~ 0.36 Hz. Since go is inversely 
proportional to the square of d, if d can be reduced to 30 nm, g 0 can be increased to 3.6 kHz. 

Supplementary Note 6: Quantum noise of optomechanical cooling 

The linearised Hamiltonian of the optomechanical system in a frame rotating at the drive 
tone frequency cud can be written as [2-6]: 

H = —hAci^a + huirJ^b — Tig(a) + a)(& + b), (7) 

where A = cud — cuo, a'(a) the creation (annihilation) operator for the cavity field variation 
and b'(b) the creation (annihilation) operator for the mechanical mode. The optomechanical 
coupling g = goVN is enhanced by the number of photons N , g 0 being the single-photon 
coupling rate. 

We obtain the Heisenberg-Langevin equations 

d{t) = [iA - a(t) + ig(b(t) + b f (t)) + 

j=e,0 

b(t) = iu m - -j-) b(t ) + ig(a(t ) + o^f)) + V7^£m(£) 
and their Hcrmitian conjugates 

d\t) = (-iA - a ] (t) - ig(b(t) + ti(t)) + y/Kj€j(t) 

j=e,0 

tf(t) = (iu) m - ) &(t) - ig(a(t) + a f (i)) + y/~fai m (t) r 

where k and y m are the decay rates of the cavity and the mechanical resonator respec¬ 
tively. The noise operators satisfy the following relationships: (£m(t)£ m (0)) = *4 <*(*), 

(Lmim = «+i)5(t) ; (£t(t)im = n e 8(t), dmim = (n e +mty (£im 0 m = 

n 0 6(t), (£o(t)£l(0)) = (n 0 + 1 )S(t)] n e K e + n 0 K 0 = n c K. 



We will rewrite this in matrix form, setting 


v(t) : = 


( a(t ) ^ 
a^(t) 
b(t) 


w(t) := 


Ej .,,0 

\ \/7m l Cm(^) / 


and 


A := 


/iA-f 

0 

ig 

\ -ig 


-.A-f 

ig 

-ig 


w 

-ig 


-iu n 


w 

-ig 

7m Q 

2 u 




*W m - \) 


thus we get 


v'(t) = Av(t) + w(t). 

Using Fourier transforms ^(/(t)) = /( u) = /j^ o /(t)e* a; * dt we have 


-iujv{ u) = Av(co) + w(co) 


which has the solution 


v(u) = (—iul — A) 1 w(u) = B 1 w(co), 


^ —i(jj 

-iA + f 

0 


0 

— 

iu> + iA 


-ig 


-ig 

V 

ig 


ig 

^1/Xc 
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-ig 

-ig ^ 

0 

1 /Xc 

ig 

ig 


-ig 

l/Xm 

0 

\ *9 

ig 

0 

l/Xm / 


-w 

-w 

ig 

ig 

ico + icu m + ^ 

0 

0 

—iu — icu m 
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For cooling we use A = — fi m and let 5 = u — D m . Applying the rotating wave 
approximation, omitting contribution from Xc, Xm, we get 

aM =-7—^- 8 

1 + 9 XcX m 

£, . Xm^/7m^m "h WXraX.c , , 

6 w =-7~r 2 -—■ 9 

1 + 9 XcX m 

From the input-output theory the output held can be expressed as 


hout Se \/ 


= 1 


Xc^e \ £ Xey^e^O £ ^y^eTmfl'XcXm £ 
1 + 9 2 XcXm) ^° 1 + 9 2 XcXm^° 1 + fi^XcXm 


Since a cancelling tone is added, the detected held is modified as 


hout Se y^e® 


Xc^e 2 Xcy^e^O 2 ^y^eTm^XcXm 2 

i + g 2 x c x™ i + g 2 XcX™ i + g 2 x c x m 




/c=e,0,m 


The spectrum analyser detects the symmetric power spectral density (PSD) [7] 
— I f / n u m a f+’mt on* 


hca 2 


e*“ (a; ut (0)a out (f) + a out (t)a^ ut (0))df. 


= 27r2^° ut ^ -a; ^ out ^ + a °«t(w)al ut (-a;)). 

Strictly speaking this is for the lab frame, however in this case only the detunings enter the 
equations, therefore we could directly substitute in Supplementary Equation 8 and 11 and 


-= Y 1 \ f ^)\ 2 C <+\\ 

k=e, 0,m ' ' 


n\ 0 being equivalent to n e>0 . Including the added noise from the amplification chain, we get 


g Ke 7m 


hLU \g2 + (f - is) - iS) 


( i 1\ |yf - i5| 2 K e K / 1\ 

f r m+ 2j + i^+(t-tf)(y-«)p r + 2) +um 


Note that the noise operators in the frequency domain satisfy the following relations: 
(H{u)£ k (u')) = 2t Tn\5{u + a;'), (ik{u)H{u')) = 2vr(n 1 fc + l)5(a; + a/). 
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The final occupation can be found via equipartition [4, 5]: 


1 + 2n m 


<^> = f 00 dc u S xx (u) 
x lp{ j- oo 2vr x 2 zp{ 

1 r 00 He ; Z" 00 1 

/ £ -mm+mme iut dt, 

"^zpf J — oo J — oo ^ 


( 15 ) 


where £(£) = x zp f(b(t) + a l (t)), x zp { being the zero-point fluctuation. Substituting in Supple¬ 
mentary Equation 9, using formulae of contour integrals and considering n 2 3> 4g 2 , /C 7 m , 7 m , 


we find 


_ «7m i . V 

U m ~ 71-n m + -- 

4^ 2 + K7 m 4g 2 + K7 n 


-n r 


1 ; C 

cT T" m + cTT" c ' 


(16) 


Supplementary Note 7: Noise from the signal generator 

For the measurement shown in Fig. 4 of the main text, we use a Phase Matrix QuickSyn 
FSW-0020 microwave signal generator. We measure a sideband noise S^ of —130 dBc/Hz 
at 120 kHz offset for a 5.1 GHz carrier signal. Cavity noise occupancy n c contributed by 
the signal’s own noise can be estimated by 


n r = 


Cm b'o * 

ftca (k/2) 2 


For a critically coupled cavity, we have 


n r = 


P in x 10 13 /Hz 
fiui 


(17) 


(18) 


The green line in Fig. 4c in the main text is plotted with the above equation. 

Measurements were also performed using an a Agilent PSG-UNY low phase noise 
option microwave signal generator which was available to us for a short time. Although 
the phase noise of the PSG-UNY at an offset of 120 kHz is specified to be -137 dBc/Hz for 
5 GHz carrier signals, the total sideband noise (amplitude and phase) was observed to be 
-132 dBc/Hz. Although we did not have time to perform a full thermal calibration of the 
setup with the PSG-UNY, we observed that the cavity noise n c was about 2 dB lower with 
the PSG generator with no evidence of additional mode heating, implying that the final 
occupation with the PSG generator would be 4.3 phonons. 
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